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ABSTRACT 

Integral  averages  of  weak  subsolutions  (and  supersolutions)  in  Rn  of 
quasilinear  elliptic  and  parabolic  equations  are  investigated.  The  important 
feature  is  that  these  integral  averages  are  defined  in  terms  of  measures  that 
reflect  interesting  geometric  phenomena.  Harnack  type  inequalities  are 
established  in  terms  of  these  integral  averages. 


AMS  (MOS)  Subject  Classifications:  35J15,  35K65 
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This  report  investigates  some  regularity  properties  of  weak  (non- 
classical)  subsolutions  of  a  wide  class  of  quasilinear  elliptic  and  parabolic 
equations  with  discontinuous  coefficients.  The  main  thrust  of  the 
investigation  is  focused  on  the  classical  concept  of  Lebesgue  point  and  the 
related  integral  average  of  weak  subsolutions.  However,  we  deviate  from  the 
classical  context  by  employing  integral  averages  that  are  defined  in  terms  of 
a  general  class  of  measures.  As  a  particular  consequence  of  the  analysis, 
Harnack  type  inequalities  for  weak  subsolutions  of  elliptic  equations  are 
obtained  that  estimate  the  difference  between  the  supremum  of  a  weak  solution 
in  a  solid  ball  in  Rn  and  its  integral  average  over  a  sphere  of  suitable 
lower  dimension,  in  the  case  of  parabolic  equations,  it  is  shown  that  the 
supremum  of  a  weak  subsolution  in  a  cylinder  in  Rn+1  can  be  estimated  in 
terms  of  its  integral  average  over  the  lateral  boundary  of  a  suitable 
subcylinder. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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Mean  Values  of  Subsolutions  of  Elliptic 
and  Parabolic  Equations 

William  P.  Ziemer 

1.  Introductory.  In  the  theory  of  regularity  of  weak  solu 
tions  that  arises  in  partial  differential  equations  and  the  cal¬ 
culus  of  variations,  the  analysis  of  the  set  of  Lebesgue  points 
frequently  plays  an  important  role.  A  point  xQ  €  Rn  is  said 
to  be  a  Lebesgue  point  of  a  function  u  if  there  is  a  number 
*,  =  Mxq)  such  that 

lim  r~n  (  | u (y )  —  i(xQ)jdy  =  0  . 

r+o  'B(xQ,r) 

Here  B(xQ,r)  denotes  the  open  ball  of  radius  r  centered  at 
Xq  .  If  u  is  a  weak  solution  of  a  partial  differential  equa¬ 
tion,  certain  regularity  properties  of  u  often  hold  at  its 
Lebesgue  points.  In  many  instances  it  is  necessary  to  investi¬ 
gate  the  Lebesgue  set  and  the  associated  integral  averages  of 
weak  subsolutions  and  supersolutions.  Perhaps  the  best  known 
results  in  this  connection  are  those  that  pertain  to  subharmonic 
functions.  These  results  state  that  if  u  is  subharmonic  in 
Rn  ,  then  u  has  the  following  sub  mean-value  properties: 

-1 

U<xn)  S  (a  (n) rn)  u  (y)  dy  , 

U  JB(x0,r) 

n  -1  t 

(1)  u(xQ)  *  <na(n)r  )  j  u  (y)  da  (y) 

9B(xq ,r) 

where  a(n)  denotes  the  volume  of  the  unit  ball  in  Rn  .  Other 
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results  in  this  direction  are  the  well  known  Harnack  inequalities 
for  weak  subsolutions  and  supersolutions  of  elliptic  and  parabolic 
equations,  [M] ,  [Tl] ,  [T2] . 

In  this  paper  we  investigate  integral  averages  of  weak  sub¬ 
solutions  of  elliptic  and  parabolic  quasilinear  equations.  This 
work  was  motivated  by  [21]  in  which  the  analysis  of  Lebesgue 
points  of  weak  subsolutions  of  parabolic  equations  played  a  vital 
role  in  determining  the  regularity  properties  of  weak  solutions 
of  a  wide  class  of  degenerate  parabolic  equations. 

In  §4  of  this  paper  we  consider  weak  subsolutions  of  equa¬ 
tions  of  the  form 

(2)  div  A(x,u,ux)  =  B(x,u,ux) 

where  A  and  B  are  measurable  functions  subject  to  certain 
structural  inequalities,  see  (12)  below.  Because  of  this  struc¬ 
ture,  it  is  natural  to  require  that  weak  subsolutions  of  (2)  lie 
in  the  Sobolev  space  W^,J)  .  Our  results  for  weak  subsolutions 
u  of  (2)  are  somewhat  analogous  to  (1)  in  that  we  relate  sup  u 
in  B(Xg,r)  to  the  integral  average  of  u  over  an  (n-k) -sphere 
of  radius  r  and  center  ,  where  1  <  k  <  p  .  This  is  a 

special  case  of  the  main  result  which  allows  consideration  of 
integral  averages  of  subsolutions  that  are  defined  in  terms  of 
certain  non-negative  measures,  (Theorem  4.6).  Of  course,  similar 
results  are  valid  for  weak  supersolutions.  In  §3,  we  consider 
u  £  Wx'p  without  assuming  that  u  is  a  subsolution  of  (2)  and 
show  that  the  integral  averages,  that  are  defined  relative  to 
certain  measures,  converge  everywhere  to  u  except  possibly  on 
a  set  whose  dimension  is  n  -  p  . 
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In  55  we  consider  weak  subsolutions  of  parabolic  equations 
of  the  form 

(3)  ufc  =  div  A{x,t,u,ux)  +  B(x,t,u,ux)  . 

We  establish  results  which  are  analogous  to  those  obtained  in  the 
elliptic  case.  In  particular  we  show  that  the  supremum  of  a  weak 
subsolution  of  (3)  in  a  cylinder  of  radius  r  in  Rn+^  can  be 
estimated  by  its  integral  average  over  the  lateral  boundary  of  a 
suitable  subcylinder  plus  an  error  term  that  tends  to  0  as 
r  ■*  0  .  Other  results  of  a  similar  nature  can  be  obtained  in  which 
the  lateral  boundary  of  the  subcylinder  can  be  replaced  by  other 
geometric  objects  of  dimension  n  . 

The  author  would  like  to  thank  Michael  Crandall  for  several 
interesting  and  helpful  discussions  during  the  preparation  of 
this  paper. 

2 .  Notation  and  Preliminaries. 

Points  in  Euclidean  n-space,  Rn  ,  will  be  denoted  by 

x  =  (x^,  x2,...,xn)  .  In  the  case  of  a  context  for  parabolic 

equations,  points  in  Rn+1  will  generally  be  denoted  by 

z  =  (x,t)  where  x  c  Rn  and  t  e  R^  .  We  will  denote  k-dimen- 

k  1 

sional  Hausdorff  measure  by  H  .  Thus,  H  is  linear  measure 

and  Hn  defined  on  subsets  of  Rn  is  Lebesgue  measure,  c.f.  [F] . 

The  Lebesgue  measure  of  a  set  E  c  Rn  is  denoted  by  !  E  J  and 
f 

!  u(y)dy(y)  will  stand  for  the  integral  average 
'E 

v  (E)'|  u(y)du(y)  ,  where  y  is  a  non-negative  Radon  measure  with 
> e  . 

u  (E)  *=  0  .  For  r:  an  open  subset  of  Rn  ,  Wi,p  ,  »>  >p  >  1  ,  will 
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denote  the  Sobolev  space  of  functions  whose  distrioutional  first 
derivatives  are  functions  that  belong  to  the  Lebesgue  space 
LPUj)  .  Wq/P($])  will  denote  the  closure  in  the  Sobolev  norm  of 

smooth  functions  with  supports  contained  in  a  .  Througnout, 
p*  =  np/n-p  will  denote  the  Sobolev  exponent  and  p'  =  p/p-1  , 
for  1  <  p  <  n  .  following  common  practice,  the  letter  C  will 
denote  a  constant  that  may  change  from  line  to  line  in  the  same 
proof. 


3.  Mean  Values  of  Sobolev  Functions, 


Whenever  E  <=  Rn  we  define,  for  1  <  p  <  n  , 


Y  p  (E)  =  inf  { 


I  V  U|  P} 


where  the  infimum  is  taken  over  all  non-negative  functions 
u  e  W1,p(Rn)  for  which  E  c  interior  {x  :  u(x)  2  1}  .  In  tne 
case  p  >  n  ,  the  definition  must  be  modified  to  require  that 
the  functions  u  e  W^'p(Rn)  have  supports  contained  in  some 
fixed  compact  set.  It  is  known,  c.f.,  [F Z ] ,  that  if  1  <  p  <  n  , 
then 

Y p (E)  =  0  implies  Hn  P+G(E)  =  0  for  each  •_  >  0 

and 

(4)  Hn  P(E)  <■  ^  implies  Yp(E)  =  0  • 
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In  the  event  that  p  =  1  ,  then  y^tE)  =  0  if  and  only  if 
Hn-1(E)  *  0  ,  [FLJ. 

3 . 1  Defuij^ion .  The  Lebes?ue  set  of  degree  s  for  a 
Lebesgue  measurable  function  u  is  the  set  of  all  points  x  e  Rn 
such  that 


lim  r-n  [  |u (y)  -  u (x)  |  s  dy  =  0  . 

r-K)  JB(x,r) 

Here  B(x,r)  denotes  the  open  ball  of  radius  r  and  center  x  . 
The  exceptional  set  of  degree  s  for  u  is  the  complement  in 
Rn  of  the  Lebesgue  set  of  degree  s  for  u  . 

It  follows  from  classical  differentiation  theory  that  in 
case  u  c  Ls(Rn)  ,  the  exceptional  set  of  degree  s  for  u  has 
Lebesgue  measure  0  .  However,  in  case  u  e  w1,p(Rn)  the  ex¬ 
ceptional  set  is  considerably  smaller.  Indeed,  we  have  the  fol¬ 
lowing,  vide  [FZ ; $ 9] : 

3.2  Theorem.  For  every  u  e  W1,p(Rn)  ,  1  s  p  <  n  ,  there 
exists  a  function  v  «  W1,p(Rn)  such  that  u  =  v  a.e.  in  Rn 
and  the  exceptional  set  of  degree  p*  for  v  has  y  capacity 
0  ,  where  p*  -  np/n-p  . 

In  particular,  this  results  implies  that  the  integral  aver¬ 
ages  of  u  converge  at  all  points  of  Rn  except  perhaps  those 
that  belong  to  a  set  of  yp  capacity  0  ;  that  is 

(3)  lim  }  u(y)dy  »  u(x) 

r;0  'B (x,r) 


-S' 
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for  all  x  £  Rn  -  E  ,  where  Yp(E)  *  0  .  The  purpose  of  this 
section  is  to  show  that  the  integral  average  taken  with  respect 
to  Lebesgue  measure  in  (5)  can  be  replaced  by  an  average  relative 
to  a  more  general  measure. 

For  this  purpose  we  introduce  the  space  of  Bessel  potentials 
g^*f  f  f  €  Lp(Rn)  ,  p  >  1  where  the  Fourier  transform  of  g^^ 
is  given  by  g^(x)  =  (l+4ir  |x|  )  .  The  space  of  Bessel  poten¬ 

tials  gx*f  ,  f  e  Lp(Rn)  is  identical  to  W1,p(Rn)  ,1  <  p  <  *>  . 
The  Bessel  kernel  is  comparable  to  the  Riesz  kernel 

R1(x)  =  |xj1_n  in  a  neighborhood  of  |x|  =  0  and  has  exponen- 

1—n 

tial  decrease  as  |x|  -*■  «  .  Moreover,  g-^  (x)  s  c|x|  for 
some  constant  C  . 

3.3  Definition.  For  each  real  number  1  ^  a  <  n  ,  let  M* 
denote  the  space  of  non-negative  Radon  measures  m  on  Rn  with 
compact  support  which  for  some  constant  M  ,  satisfy  the  growth 
condition 


m(Br)  =s  Mr0"0 

on  open  n-balls  of  radius  r  . 

If  p  >  1  and  m  is  a  non-negative  Radon  measure  such  that 

]|  *zn  ||  ,  <  »  ,  where  p'  =  p/p-1  ,  then  clearly 

x  Lp  (Rn) 

m  €  [W1,p(Rn)]*  ,  the  dual  of  W1,p(Rn)  .  To  see  this,  let 
u  <  W^'p(Rn)  be  written  as  u  *  g^*f  ,  f  £  Lp(Rn)  .  Then,  by 
changing  the  order  of  integration,  we  have 
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and  thus,  }  udm  defines  a  bounded  linear  form  on  W  (R  )  • 
Again,  by  changing  the  order  of  integration,  we  have 

J  (g^m)*'  dx  *  |  g1*(g1*m)1/p"1dm  . 

Therefore,  by  appealing  to  Theorems  3.3  and  3.5  of  [AMI,  it  fol-^ 

lows  that  if  p  >  1  ,  1  *  <*  <  P  and  m  6  Ma  '  then  *  *  >} 

In  case  p  -  1  ,  it  follows  from  Theorem  4.7  of  IMZ]  that 
m  £  lW1'1(Rn)l*  if  m  e  M+  .  The  main  results  of  [MJ  and 

[MZl  thus  yield  the  following. 


3.4  Theorem.  Let  «  c  Rn  be  a  bounded  Lipschitz  domain 
and  suppose  1  s  P  <  n  .  Let  m  be  a  non-negative  Radom  measure 

supported  on  0  with  m  e  M+  if  P  -  1  223.  m  €  K  if  1  S  ° 
Th»r»  is  a  constant  C  =  C(0,n,p,M)  such  thatj-f  u  c  WlfPlfl)  * 

then 

,{  |u(y)  -  {  udm|p*ay)1/p*  *  C(ja  |vu|p>1/p  . 

Thl.  result  will  be  crucial  in  e.tabliehing  (5)  where  Lebesgue 
measure  in  the  integral  average  is  replaced  by  »  •  • 

If  u  €  W1,p(Rn)  .  define  a  measure  yu  by 


«L 


Yu(E)  =  |7u|Pdx  . 

E 

whenever  E  is  a  measurable  subset  of  Rn  . 
Let 


P  *  Rn  n  (x  :  lim  sup  rP”n  y  [B(x,r)]  >  0}  . 
p  r  +  0  u 

3.5  Lemma.  I£  ucW^,p(Rn)  ,  1  s  p  <n  ,  then 
-  0  • 

Proof :  Let  *  Rn  n  (x  :  lim^gup  rp  n  yu[B(x,r]  >  k-1 } 
for  each  positive  integer  k  .  It  follows  from  [F; §2.10.19] 
that  there  is  a  constant  C  such  that 


H"'P<V  S  <*  T^ljy  . 


Note  that  yu  is  absolutely  continuous  with  respect  to  Lebesgue 

measure.  Now  y^A^.]  <  »  ,  Hn_p(AJt)<  »  ,  Hn(AJ<)  =  0  , 

vu(Ak*  =  0  and  Hn-1(Ak)  =  0  .  Hence  Hn”P(pn>p)  =  0 

because  P  =  uA,  . 

n-p  Tc 

Let  x0  *  pn-p  and  consider  B(xQ,l)  .  If  u  e  W1,p(Rn)  , 
1  s  P  <  n  ,  and  if  m  is  a  measure  satifying  the  hypothese  of 
Theorem  3.4  where  we  set  ft  =  B(xQ,l)  ,  then 


(6)  (|  |  u  (y ) 

JB(xq,1) 


udmjp*  dy)  */p*  * 


| 7U |P) 1/P 

l) 
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Let  Tr  :  B(x0,l)  ■>  B(xQ,r)  be  defined  by  Tr(xQ+y)  =  xQ  +  ry 
where  y  e  B(0,1)  .  If  we  define  ur  =  u  o Tr  ,  then  (6)  implies 
that 


f 

(  |u  (y)  -  u  •  T  dm 

Jb(xq,1)  r  J  r 


dy) 


1/P* 


s  C( 


7u 


B(x0,l) 


Pj  1/P 


which  is  the  same  as 


(7) 


<r'n  } 


B(xQ,r) 


u(y)  ~ 


u  o  Tr  dm  |P*  dy) 1//p 


s  C(rp“n  [  |7u|p)1/p  . 

JB(xA,r) 


Because  x„  P  _  ,  the  right  side  of  (7)  tends  to  0  as  r  4-  0 
*  u  n-p 

From  Theorem  2.2,  there  is  a  set  E  with  Yp(E)  =  0  such  that 

(8)  lim  r“n  f  |u(y)  -  u(xn)|p*dy  =  0 

r+0  JB(xQ,r) 


whenever  xQ  fc  E  .  Note  that  (4)  and  Lemma  3.5  imply  y  (P  )  =0 
Now  by  setting  A  =  E  u  Pn_p  »  the  following  result  follows  from 
(7)  and  (8). 


3.6  Theorem.  Suppose  u  t  W1,p(Rn)  ,  1  <.  p  <  n  .  There  is 
a  set  A  c  Rn  with  Yp(A)  =  0  such  that  if  xQ  4  A  and  m  is 
a  non-negative  Radon  measure  supported  on  B(xQ,l)  with  m  e  M+ 
if  p  *  1  and  m  e  irt*  if  1  s  a  <  p  ,  then 


(9) 


lim 
r  +  0 


u  o  Tr  dm  =  u  (Xg)  . 
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Of  course,  the  most  interesting  measures  to  consider  are 
those  that  passess  some  homogeneity  properties.  For  example, 
if  m  =  H11”1  I  3B  (Xg ,  1)  ,  then  clearly  m  c  and 


u  o  Tr  dm  =  /  u (y)dHn_1 

*  3B (xrt ,r ) 


Thus,  we  have 


3.7  Corcdlar^.  If  u  e  W^'^{Rn)  there  is  a  set  A  with 
(A)  =  Hn-1(A)  =  0  such  that 


lim  i  u(y)dHn-1(y)  =  u  (x  ) 

r+0  J  3B(xn,r) 


whenever  xQ  A  . 

A  similar  result  could  be  obtained  by  taking 
m  =  Hn  1  7?n  ^(Xg)  where  irn  *(Xg)  is  an  (n-1) -dimensional 

plane  passing  through  Xg  .  In  the  event  that  u  e  W^'^(Rn)  , 

1  <  p  <  n  ,  let  m  =  Hn  R  Sn  k(xQ,l)  where  1  £  k  <  p  is  an 
integer  and  Sn  k(Xg,l)  is  an  Sn-ksphere  with  radius  1  centered 
at  Xg  .  Then  we  have 


lim 

r*° 


u(y)dHn"k{y) 


U(Xg) 


for  y  -  q.e.  xQ  e  R  .  As  in  the  case  p  =  1  ,  a  similar  re¬ 
sult  follows  by  taking  m  =  H*1”14  ffn"’k(Xg)  . 
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4 .  Lebesgue  points  of  subsolutions  of  elliptic  equations 

In  this  section  we  obtain  results  for  weak  subsolutions  and 
supersolutions  of  quasilinear  elliptic  equations  with  measurable 
coefficients  that  are  analogous  to  Theorem  3.6  and  its  corollaries. 
In  the  case  of  weak  subsolutions  we  will  show  that  (9)  holds  at 
all  points  Xq  and  more  importantly,  that  the  sup  u  in  B(Xg,r) 
is  bounded  above  by  the  left  hand  side  of  (9)  plus  a  term  that 
tends  to  0  as  r  +  0  . 

The  equations  that  are  considered  in  this  section  are  of  the 

form 

(11)  div  A{x,u,ux)  =  B(x,u,ux) 

where  A  and  B  are,  respectively,  vector  and  scalar  valued 

1  . 

Baire  functions  defined  on  a  x  R  x  R  .  Here  n  is  an  open 
subset  of  Rn  .  The  function  A  and  B  are  required  to  satisfy 
the  following  structural  inequalities: 

|A(x,u,w)  |  s  a0|w|p“1  +  a^Ju^-1  +  a2 

(12)  |  B(x,u,w)  |  s  t>0jw|p  +  bjJ  w|  P”^  +  b2|u|p-1  +  b3 

A(x,u,w)  •  w  2:  |w|p  -  CjJujP  -  C2  . 

We  assume  1  <  p  <  n  ,  and  Sq  ,  bg  are  non-negative  constants. 

The  results  below  are  valid  if  the  remaining  coefficients  are 
non-negative  measurable  functions  that  are  assumed  to  lie  in 
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appropriate  L ^(q)  spaces,  vide  [GZ]  ,  but  to  minimize  technical 
detail,  we  will  assume  that  the  coefficients  are  bounded  by  some 
constant  K  >  0  : 


(13)  ai(x)  *  K  ,  b^x)  *  K  ,  c^x)  ^  K  . 

A  function  u  e  wioc^  is  calle<*  a  weak  subsolution  (supersolu- 
tion)  of  (11)  if 

(14)  I  A(x,u,ux)  •  V*  +  B(x,u,ux)^  £  0  (2  0) 

ft 

for  all  bounded  <t>  *  0  ,  <t>  e  Wq'p(^)  . 

The  following  result  is  due  to  Trudinger  [Tl]  whose  proof  is 
based  on  the  familiar  Moser  iteration  method,  [M] .  Therefore,  a 
complete  proof  will  not  be  given;  however,  for  the  convenience 
of  the  reader,  we  will  give  an  outline  and  provide  the  main  steps 
of  the  proof. 

Let  u  e  be  a  weak  subsolution  of  (11)  that  is 

bounded  above  by  L  in  ft  .  Choose  Xg  e  ft  and  for  each  real 
number  k  ,  let  =  (u-k)+  and  define 

wk(r)  =  sup  (uR<x)  :  x  e  B(xQ,r)}  , 

(15) 

y(r)  =  sup  (u (x)  ;  x  c  B(xQ,r)}  . 

4.1  Theorem.  Let  u  e  ^ioc^  ,  1  <  p  <  n  ,  be  a  weak  sub¬ 
solution  of  (11)  such  that  0  s  u(x)  s  L  for  each  x  e  c  .  There 
is  a  constant  C  depending  only  on  L  ,  n  ,  p  and  the  structure 


2- 


(12)  such  that  if  B(xQ,r)  c  u  ,  then 

(J-  |  u  (r)  -  u|y)1/y  s  C[v  (r)  -  y(  r/2)  +  a(r)] 

JB(x0,7r/8) 

where  a(r)  *  r  +  Kr  +  (Kr)p/p_1  ,  and  y  <  n(p-l)/(n-p)  . 

Proof :  Let 

-bnu 

(16)  <j>  =  nPe  (y(r)  +  a(r)  -  u) 

where  n  €  Cg[B(x0,r)]  and  a  <  0  .  Then  4,  is  a  bounded,  non¬ 
negative  test  function  that  can  be  employed  in  (14) .  If  we  set 

v  *  u (r)  +  a(r)  -  u 

the  structure  (12)  and  elementary  estimates  yield 

(17)  |  nPva-1|vv|p  s  C (a )  |  (nP  +  |Vn|P)vp+a+1 

where  C(a)  is  a  finite  constant  when  a  is  bounded  away  from 
0  and  1  -  p  .  Let  w  *  vq  where  pq  »  p  +  a  -  1  and  apply 
Sobolev’s  inequality  to  find  that  (17)  yields 

(18)  llnw^p  s  C (a) jj  (n  +  |Vn|)w|^ 

where  0  *  n/n-p  .  For  0  <  s  <  t  s  r  ,  let  n  €  CptBCx^rt)]  be 
such  that  0  s  n  s  1  ,  n  =  1  on  B(xQ,s)  and  |Vn|s  2(t~s)  1  . 
Then  (18)  implies 


U9>  l|w||op),  s  C(»)(t-s)"p||w||p.t 

Let  Tj  »  r(2~*  +  2-^-  j  for  j  *  0,1,... ,  and  iterate  inequality 

(19)  to  obtain  for  any  pQ  >  0 

(20)  e(-».r/2)  ^  Ce(-p0,3r/4) 

where 


0(P,r)  =  (P  vp)1/p  . 

J  B(xQ,r) 

By  performing  a  finite  iteration  of  (19)  and  choosing  the  r^  in 
a  different  way,  we  have 

(21)  e (Y,7r/8)  s  C8(p0,3r/4) 

for  any  pQ  >  0  and  Y  <  n(p-l)/(n-p)  .  Return  now  to  the  de¬ 
finition  of  the  test  function  4  and  notice  that  in  case  a  ■  -1  , 
substitution  of  4  into  (14)  gives 

J  np|  v  logv|  p  s  C  |  |  Vn  |  p 

whenever  n  e  CQ[B(xQ,r)]  .  Thus,  from  the  John-Nirenberg  lemma, 
c.f.  [GT,  p.158],  it  follows  that  there  exist  constants  C  and 
Pq  >  0  such  that 


B(xq,s) 


-po  f 

v  J 


B(x0,s) 


p0 

v  w  <;  CS 


2n 
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for  0  <  a  s  3r/4  .  That  is, 


(22)  0(po,3r/4)  *  e(-p0,3r/4)  . 

Thus,  (20) ,  (21)  and  (22)  imply 

0(Y,7r/8)  £  C0(-«,r/2) 


or 


({-  lu(r)  +  a (r)  -  up)1/*  *  (f-  v*)1/y 

'B(xQ,7r/8)  ’ b(xq ,7r/8) 

s  C  min  v  «  C(p(r)  -  n(r/2)  +  a(r)]  . 

B(x0,r/2) 

4.2  Remark.  Note  that  the  right  hand  side  of  the  inequality 
in  Theorem  4.1  tends  to  0  with  r  and  therefore 

lim  /  u(y)dy  *  lim  y(r)  *  lim  sup  u(x)  . 

r+0  ;B(x0,r)  r+0  x-*-x0 

This  implies  that  every  bounded  weak  subsolution  of  (11)  is  upper 
semicontinuous  on  a  after  redefinition  on  a  set  of  measure  0  . 

If  bg  »  0  in  the  structure  (12) ,  then  the  same  conclusion 
will  hold  without  assuming  that  u  is  bounded  on  a  .  In  this 
case  it  is  possible  to  show  that  a  weak  subsolution  u  e 
is  locally  bounded  above  on  n  ,  [SEJ,  On  the  other  hand,  u  is 
not  necessarily  bounded  below  and  therefore,  the  estimate  in 
Theorem  3.1  may  not  hold.  However,  one  can  easily  show  that  for 
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9 


U 


each  k 


R 


1 


lim  1 

r+0  >  B (xQ , 7r/8) 


|uk<r>  -ukiv 


0  . 


To  see  this,  replace  u  by  and  y(r)  by  yk(r)  in  the 

definition  of  <f>  ,  (16)  ,  and  observe  that  the  proof  of  Theorem 
3.1  yields 


(23)  (/  Iliv(r)  -u.  |Y)1/Ysc(k)  [y.  (r)  -  y,  (r/2)  +a(r)]  . 

' B (xQ , 7r/8)  K  K  k  k 


Let  A(xn)  *  lim  y(r)  and  choose  k  <  A (x.)  .  Then 
u  riO  0 


y,  (r)  -us  ( A  (x- )  -  k)  |  B  (x_  ,  7r/8)  n  {u  <  k}  | 
B(x0,7r/8)  k  k  0  0 


and  thus,  it  follows  from  (23)  that 


lim  r  n|s(xn  ,  r)  n  (u  <  k}|  =  0 
r-K)  u 


Therefore,  for  each  t  >  0 


lim  r  n|B(xn  ,  r)  n  {x  :  |u(x)  -  A (xft) |  >  e } |  =0  . 
r+0  u  u 

This  states  that  u  is  approximately  continuous  at  xQ  and  that 
its  approximate  limit  at  Xq  is  A(Xq)  •  However,  every  measur¬ 
able  function  is  approximately  continuous  almost  everywhere  and 
therefore,  every  weak  subsolution  of  (11)  with  bQ  =  0  in  (12) 
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«* 


is  upper  semi continuous  on  ft  after  redefinition  on  a  set  of 
measure  0. 

We  will  now  show  that  (9)  holds  at  all  Xg  e  ft  if  u  is  a 
subsolution  of  (11). 

4.3  Theroem .  Let  u  e  ,  1  <  p  <  n  ,  be  a  weak  sub¬ 

solution  of  (11)  such  that  0  <,  u(x)  s  L  for  each  x  €  ft  .  For 
each  xQ  e  ft  with  B(Xg,rg)  c  ft  for  some  rQ  >  0  ,  let  m  e  M*  , 
1  z  a  <  p  ,  be  a  Radon  measure  supported  on  B(Xg,rg/2)  .  There 
is  a  constant  C  depending  only  on  M  ,  L  ,  n  ,  p  and  the  struc¬ 
ture  (12)  such  that 


(24) 


(f  |u(y)  -  [  u  o  T  dm|p*dy)1</p* 

JB(x0,r)  j  r 

s  C [y (r)  -  y (r/2)  +  a(r))p"1/p 


for  each  r  <  rQ/2  where  Tr  :  B(Xg,rg/2)  •*  B(xQ,r/2)  is 
defined  by  Tr(xQ+y)  ■  xQ  +  r/rQy  . 


Proof .  We  proceed  to  obtain  an  estimate  of 


|  vu| 


2 


JB(xQr) 

as  in  [GZ ] .  Let  n  €  Cg [B(Xg,7r/8) J  be  a  cut-off  function  such 

bgii 

that  n  =  1  on  B(Xg,r/2)  and  let  $  *  np  e  u  .  Hence,  by 
substituting  <j>  into  (14)  and  utilizing  the  structure  (12)  ,  it 
follows  that  there  is  a  constant  C  such  that 
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(25) 


nP | vu | p  s  c[j  nP~1 1 v n | P_1 1 Vn |  +  j  nP 

+  |  np  ^ | V n |  +  |  np|vu|p  . 

As  in  the  proof  of  Theorem  4.1,  let  v  =  p (r)  -  u  +  a(r)  and  use 
(17)  and  Theorem  4.1  to  obtain,  for  sufficiently  small  e  >  0 

(26)  np_1i  vu  | P_1  f  vn  |  =  j  ( np_1v“  (1+e)  /p'  |  7u  I15”1)  (v(1+e)  /p'  |  Vn  | ) 

*  ({  npv‘(1+e,|w|P)P-1/P  (J  v(P-1)(1+e)|7n|p)1/p 

^  C(  (nP  +  |7n|P)vp“1‘e)p“1/p 

([p(r)  -  p (r/2)  +  a(r)] (p"1) (1+e)  •  rn“p) 1/p 
s  C([p(r)  -  p  (r/2)  +  a(r)]  (P"1"e)rn_P)p";L/P 
( [u (r)  -  p (r/2)  -  a(r)] (P”1) (1+e)rn"p)1/p 
s  C [p (r)  -  p (r/2)  +  a(r)]p_1rn"P  . 

Also,  the  last  term  in  (25)  can  be  written  as 

-  |  (n|vu|)P_1„ 

s  %  j  nP|»ulp  +  c  |  . 

Hence,  it  follows  from  (25)  and  (26)  that 
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t. 


D-n 


vu 


B (xQ ,r/2) 


|p  s  Ctw(r)  -  u  (r/2)  +  atr)]15”1 


Now  apply  Theorem  3.4  as  in  the  proof  of  (7)  to  establish  the  con¬ 
clusion. 

Because  the  right  hand  side  of  (24)  tends  to  0  as  r  +  0  , 
the  following  is  immediate. 


4.4  Corollary.  If  u  e  W^p(ft) 
(11)  on  ft  ,  then 


is  bounded  subsolution  of 


lim 
r  +  0 


u  e  Tr  dm  *  u(Xg) 


for  each  xQ  e  ft 

4.5  Remark.  As  in  (10),  if  we  let  m  =  Hn-k  |  Sn-k (xQ,r0/2)  , 
we  have  for  each  xQ  e  ft 


lim 

r+0 


<Vr) 


u(y)dHn~k(y) 


u(x0) 


whenever  k  is  an  integer  such  that  1  s  k  <  p  .  A  similar 
result  holds  if  m  is  taken  as  the  restriction  of  H  to  an 
(n-k) -plane  passing  through  x^  . 

Finally,  we  establish  an  inequality  for  weak  subsolutions  of 
(11)  which  is  reminiscent  of  the  weak  Harnack  inequalities  proved 
in  [Tl]  . 
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be  a  weak  sub- 


% 


4.6  Theorem.  Let  u  e  W^p(fi)  *  1  <  P  <  n  , 
solution  of  (11)  and  assume  that  0  <  u  <  L  on  fi  .  For  each 
xQ  e  fi  with  B  (Xg  ,rg)  c  $2  for  some  r^  >  0  ,  there  is  a  non¬ 
negative  function  g(r)  with 


0 


g(r) 


P' 


dr 

r 


oo 


such  that 


supu  s 
B(xQ,r) 


Sn~k(xQ,r) 


u(y)dHn“k(y)  +  g (r) 


for  0  <  r  <  r  /2  .  Here  k  is  an  integer  such  that  1  <  k  <  p 


Proof.  This  follows  immediately  from  Theorem  4.1  and  (24) 
for,  setting  c(r)  = 


Sn"k(xQr) 


u(y)dHn  k(y)  ,  we  have 


y(r)  -  c (r) 


(xA# 


7r/8) 


u(r)  -  c (r) 


•I 


B(xQ,7r/8) 


v  (r) 


-  u|  +  | 


B(xQ,7r/8) 


u  -  c (r) 


s  C[y(r)  -  u(r/2)  +  a(r)]p_1/p  . 


Now  set  g(r)  =  C[u(r)  -  w(r/2)  +  a(r)]p  ^p  and  the  conclusion 
follows . 


5.  Parabolic  Equations.  In  this  section  we  consider  parabolic 
equations  of  the  form 
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(27) 


ufc  *  div  A(x,t,u,ux)  +  B(x,t,u,ux) 


where  A  and  B  are  Baire  functions  defined  on  *  R1  *  Rn  . 
Here  ■  n  *  (0,T)  where  fi  is  an  open  subset  of  Rn  . 

The  structure  imposed  on  A  and  B  is  similar  to  that  con¬ 
sidered  in  the  elliptic  case#  (12): 

|A(x,t,u,w)  I  i  aQ|w|  +  ajuj  +  a2 

(28)  |B(x,t,u,w)|  s  bQ ! w 1 2  +  bx|w|  +  b2|u|  +  b3 

A(x, t ,u,w)  •  w  2  1  w |  -  c^Juj  "  c2  . 

aQ  and  bQ  are  non-negative  constants  and  the  remaining  coef¬ 
ficients  are  required  to  lie  in  appropriate  Lebesgue  spaces,  c.f., 
[LSU] ,  [T2] ,  but  as  in  S4,  for  simplicity  of  exposition,  we  will 

assume  that  all  coefficients  are  bounded  by  some  constant  K  >  0  . 

1  2 

A  function  u  t  W*^(0T)  is  called  a  weak  subsolution  (super¬ 
solution)  of  (27)  if 

(29)  |  -  u$t  +  A  •  V*  -  B*  s  0  (20) 

for  all  bounded  <fr  2  0  ,  4  «  wJ'2(fiT)  . 

The  object  of  this  section  is  to  establish  results  for  sub¬ 
solutions  of  (27)  analogous  to  those  in  Theorems  4.3  and  4.6.  For 
this  purpose  we  consider  an  arbitrary  point  ZQ  e  nT  ,  and  for  con 
venience,  set  Zg  *  (0,0)  .  We  shall  utilize  space-time  cylinders 
of  the  follwoing  form: 
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R(r)  =  B (r)  *  (o^^r2) 


R~ (r)  =  B(ar)  x  {t  ^r:2 , 1 2r2 ) 


R*(r)  =  B(er)  x  {P;Lr2,p2r2) 


where  B(r)  denotes  the  n-ball  of  radius  r  and  center  =  0 
and  where 


0<  8  <  a  <  1  , 


(30) 


<  Pl  <  p2  <  <  0  <  p4  <  o2 


Note  that 

(31)  R(r/2N)  c  R~  (r) 

1  2 

for  some  integer  N  .  Finally,  for  u  e  W^c(nT)  a  subsoluti°n 
of  (27) ,  let 


y(r)  =  sup{u(z)  !  z  e  R(r) }  . 


The  following  result  is  proved  in  [T2]  and  its  proof  runs  parallel 
to  that  of  the  elliptic  version  in  Theorem  4.1. 

1  2 

5.1  Theorem.  Let  u  «  W^c(nT)  be  a  weak  subsolution  of 
(27)  such  that  0  s  u(z)  s  L  for  each  Z  e.  0 _  . 


There  is  a 


number  y  >  1  and  a  constant  C  depending  only  on  L  ,  n  ,  and 
the  structure  (28)  and  (30)  such  that  if  R(r)  c  ,  then 

(i  |u(r)  -  u|VA  s  C  min  (u(r)  -  u  +  a(r)] 

^  R*  (r)  R~(r) 

where  a(r)  =  r  +  Kr  . 

2 

We  now  will  establish  an  estimate  on  the  growth  of  V  -norm 
of  u(r)  -  u  where  u  is  a  weak  subsolution  of  (27).  This 
estimate  is  similar  to  the  one  that  appears  in  [Z2,  Theorem  4.3] 
and  is  the  parabolic  analogue  of  the  estimate  that  appears  in 
Theorem  4.3  above. 

For  this  purpose  let  u  be  a  bounded,  non-negative  weak 
subsolution  of  (27)  and  let 

2  bnu 

(32)  $  =  n  u  e  0 

where  n  is  a  smooth  cut-off  function  that  will  be  specified 
below.  Because  u  is  assumed  to  be  bounded,  we  may  assume  that 
the  terms  involving  | u |  ,  b2|u|  ,  and  c^Ju]  in  (28)  can  be 

absorbed  in  a2  ,  b^  and  c 2  respectively.  Therefore,  substi¬ 
tution  of  (32)  into  (29)  yields 

f  2  bnu  2  bnu  7 

(33)  J  n  e  uut  +  n  e  (1  +  bQu)(j7u|  -  c2) 

^  ||  2nu  e  0  [a0)7u|  +  a2]|7n| 

+  ||  n2  e  0  u(bQ|7u|2  +  b1i7u|  +  b3J 
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Now  let  v  =  y(r)  -  u  (for  simplicity  we  will  write  y  =  u(r)  ) 
and  define 


(34) 


>  .-1  bou,,-i  .  .-1  b0y,.-l  , 

f(u)  =  bQ  e  (bQ  -  u)  -  bQ  e  (bQ  -  u)  . 


Then, 


(35) 


bou 

f ’ (u)  »  -u  e  , 


b0u 

f"(u)  =  -e  u  (1  +  bQu)  , 


^(y2  -  u2)  <  f(u)  ^  y  e  0  (y  -  u)  for  0  <  u  <  y  . 


There  is  a  constant  C  such  that  (33)  can  be  written  as 


(36) 


’[  2  9f(u)  f 

J  3t  J 


'  n2 | Vu | 2 


sC[l  n ! Vu  j  | v  n  !  + 


n2  + 


n | Vn |  +  I  I  n2 I Vu  j ] 


The  constant  C  depends  upon  the  structure  (28)  and  the  bound 
for  |u|  .  Now  define 


R* (r)  =  R*(Yr) 

where  0  <  y  <  1_  and  choose  n  so  that  r,  =  1  on  R'  (r)  and 
spt  n  c  R*(r)  .  Refer  to  (30)  and  set  t^  *  p^r2  and  t2  =  P2r2 
Choose  t*  e  (t1»t2)  so  that 
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Hi 


j  n2  (x,  t*)  f  2  (x,  t*)  dx  >  \  sup  f  r,2  (xft)  f2  (x,  t)  dx  . 
'  (tX,t2)  J 


If  we  replace  t  by  t*  in  (36),  we  are  led  to  the  following 
estimate: 


(37) 

sup 

[nf(x,t)]2dx  +  [f  n2(7v|2dxdt 

<vv 

J ; 

£  C[ 

nj  Vu)  |  7n)  +  |  n2  +  |  n  |  v  n  | 

+ 

nfc|vuj  +  |  n[nt|f]  . 

We  now  proceed  exactly  as  in  the  proof  of  Theorem  4.3  above  to 
estimate  the  rate  at  which  the  right-side  of  (37)  tends  to  0. 
Indeed,  by  employing  (31)  and  Theorem  5.1,  we  find  that  there  is 
a  constant  C  depending  only  on  the  given  data  such  that  the 
first  four  terms  on  the  right-side  of  (37)  are  bounded  by 

(38)  C  rnlv(r)  -  u(r/2N)  +  a(r)]  . 

Refer  to  (35)  to  find  that  |f(u)|  s  C^|v|  and  therefore,  by 
using  Theorem  5.1  again,  we  have 

(39)  ||  n|nt!Jf(u)|  s  C  rn(y(r)  -  »(r/2N)  +  a(r)]  . 

Referring  again  to  (35),  the  following  lemma  now  follows  from 
(37),  (38)  and  (39). 


1,2 


3.^  hec  u  e  wioc be  a  bounded ,  non-negative 

weak  subsolution  of  (27) .  Then, 


sup 

(tl't2) 


B  (3  Y  r) 


[y  (r)  -  u]  dx  + 


Mr*  (r) 


Vu  |  2  dxdt 


<  C  rn[g(r)  -  y (r/2N)  +  a(r)] 


whenever  R(r)  c  . 

We  are  now  in  a  position  to  establish  the  parabolic  analogue 
of  Theorem  4.3.  This  result  is  concerned  with  the  behavior  of  a 
weak  subsolution  in  a  neighborhood  of  an  arbitrary  point  ZQ  e  «T  . 
We  will  continue  to  assume  that  Zq  =  (0,0)  and  the  geometric 
configuration  imposed  by  (30). 


5.3  Theorem.  Let 
weak  subsolution  of  (27) 

for  some  rQ  >  0  ,  let 
supported  on  B  (3  y  )  * 
only  on  the  bound  for  u 


1  2 

u  c  wiqC(°t)  be  a  bounded,  non-negative 
.  For  each  ZQ  e  nT  with  Rz^(rQ)  c 

m£M^,lsa<2,  be  a  Radon  measure 
(t)  .  There  is  a  constant  C  depending 
and  the  given  data  such  that 


|u(x,t)  -  I  u  o  Tr(y,t)  dm(y) 


2  dxdt 


s  C  [y (r)  -  y(r/2N)  +  a(r)] 


for  each  r  s  rQ  where  Tr  :  B(gYr0)  x  { t >  -►  Bteyr)  x  {t}  is 
defined  by  Tr(x,t)  =  (r/rQx,t)  . 
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1  2 

Proof.  Because  u  e  wioC  ^T^  '  an  elementarY  fact 

1  2 

concerning  Sobolev  functions  that  u(*,t)  e  *  {t})  for 

a.e.t.  For  all  such  t  apply  Theorem  3.4  to  obtain 


B(0  yr) 


|U(X,t)  - 


u  o  Tr (y,t)  dm(y) |  dx 


^  C  r 


2-n 


B  (3  y  r) 


7u  (x / 1)  |  dx  . 


Now  integrating  with  respect  to  t  from  y^  =  yp^r  to 
2 

yt2  =  yp2r  and  applying  Lemma  5.2  yields  the  desired  result. 


5.4  Corollary.  For  each  ZQ  e  Si T  with  Rz  (r^)  c  nT  for 
some  r.  >  0  ,  there  is  a  non-negative  function  g(r)  with 


g(r) 2  ^  <  » 


such  that 


,yt2  , 

supu  s  /  u  «  T 

R(r)  Jyt,  *  r 


(y,t)  dm(y)  dt  +  g(r) 


for  0  <  r  5  r^  . 


Proof.  From  Theorems  5.1  and  5.3,  it  follows  that 
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I  U  o  Tr(y, 


t)  dm(y) 


L-! 


u(y ,  t) dHn_1 (y) 


S  (xQ ,3  yr)  *{t} 


and  therefore 


n 

JYt, 


u  •  Tr(y,t)  dm(y)  dt 


-  f 2 1 

Jyt,  ln-1 


u(y,t)dHn  1 (y) 


Yt"l  S“  *(x0,8yr)x{t} 


3*R'  (r) 


where  3*R' (r) 


u(z)dHn{z) 

denotes  the  lateral  boundary  of  the  cylinder 


R’ (r)  . 

We  have  thus  proved 


5.5  Corollary.  For  each  e  with  Rz0(r0)  c  nT  for 
some  rQ  >  0  ,  there  is  a  non-negative  function  g(r)  with 

f  g(r)2  ^  «  - 
Jo  r 

such  that 

supu  z  /  u(z)dHn(z)  +  g(r) 

R(r)  i  3*R’ (r) 


for  0  <  r  s  rQ 

As  in  S3  and  S4,  one  could  establish  similar  results  by 
making  different  selections  for  the  measure  m  .  For  example,  in 
addition  to  the  choice  made  in  the  above  corollary,  one  could 
also  take  m  *  Hn  1  j  nn  ^(Xq)  *  {t>  where  nn-1  is  an  (n-1)- 
plane  that  is  orthogonal  to  the  t-axis. 
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